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SUMMARY 

A method is presented for determining, hy stej>-'by— step integra- 
tion, the trajectories of water drops aroijnd any body in two- 
dimensional flow for which the streamline velocity components ane 
known or can he computed.. The method is general and considers the 
deviation of the water drops from Stokes' law because of speed and 
drop size. 

The equations are presented in general form and then, to 
illustrate the procedure, water-drop trajectories are calculated 
about a 12— percent— thick symmetrical Joukbwski profile chosen to 
simulate an NACA 0012 section. 

The method provides a means for the relatively rapid calculation 
of the trajectory of a single drop without the utilization of a 
differential analyzer. 

In addition, consideration is given to the ma.YiTni]Tn possible rate 
of water-drop impingement on a body. 


INTRODUCTION 

O 

As part of a comprehensive research program directed toward a 
fundamental understanding of the mechanism of thermal ice— prevention 
for airplanes, the NACA has undertaken an investigation of the 
performance of a heated wing or empennage section in specified icing 
conditions. One of the first essentials of such an investigation is 
a method, for estimating or calculating the area over which water will 
strike the wing, and the distribution of water impingement over that 
area. 
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Previous investigations of this nature have considered the 
vater-drop trajectories about cylinders (references 1 and 2) and 
those about an airfoil (reference 2). In both of these studies the 
assvunption was made that the water drops were sufficiently small to 
obey Stokes' law of resisteince. A more recent investigation for 
cylinders '(reference 3) shows that, for the velocities of airplanes 
and the drop sizes frequently pre.sent in clouds, Stokes' law does 
not apply and the force acting on the drop can be determined only 
from a knowledge of the drag coefficient for spheres. 

In reference 3> ' "therefore, the methods of previous investigators 
are modified to include more, accurate, values of the drop drag coeffi- 
cient. Equations are established for the determination of the drop 
trajectories around cylinders, spheres, and ribbons. These equations 
are derived in forms suitable for use ,.in a differential analyzer and, 
by utilizing an analyzer, the trajectories are calculated and plotted. 

The designer of a heated wing, desiring to know the rate and 
area of water iH5)ingement on the leading edge in a specified cloud at 
a given flight speed,, might assume that the in5)ingement will be 
identical to that for a cylinder with a radius equal to the wing 
leadiiig— edge radius. There is some question, however, as to the 
accuracy of this assumption for the larger drop sizes and for wing 
sections with small leading-edge radii. The designer is therefore 
confronted with the desirability of employing a computation method, 
preferably without the necessity of a differential analyzer, which 
will provide some indication of the area and distribution of water 
impingement . 

In this report, the water-drop trajectory eqxiations used in 
reference 3 have been modified to establish a step-by— step integra- 
tion method applicable to any two-dimensional flow for which the 
streamline velocity components are known or can be approximated. If 
drop trajectories for a large range of airspeeds and drop diameters 
are required, the computation time is large, and the desirability of 
access to a differential analyzer becomes evident. The integration 
method presented herein, however, does permit the calculation of 
any desired number of trajectories without resort to an analyzer; 
it also provides a means for estimating the error which will be 
incurred by replacing the airfoil by a. cylinder with a radius eqiial 
to the leading-edge radius of the airfoil. In addition, water-, 
impingement data over the entire airfoil surface can be obtained by 
the integration method. 
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SYMBOLS 

Coefficients 

radius of a drop, microns ( 3*28 x 10“® ft) 

average acceleration of a drop over an interval of time, 
feet per second per second 

acceleration of a drop at the heginning of an interveil 
of tine, feet per second per second 

dimensionless ratio \ ' 

V2^K/ 

chord of airfoil, feet 

dimensionless concentration factor 

drag coefficient for spheres 

efficiency of deposition of water drops on a body 

dimensionless drop— inertia q^uantity, defined in equation (4) 

dimensionless drop— inertia quantity for which no drops 
impinge upon a body 

unit length of span of an airfoil 

intensity of water-drop impingement at a point on the surface 
of a body, pounds per hour, sqiiare foot 

rate of water-drop impingement on a body, pounds per hour, 
foot span 

nondime ns ional distance 

local Reynolds number of a drop relative to the air stream 

Reynolds number of a drop moving with free— stream velocity 
in still air 

disteince along the surface of a body from zero percent 
chord, feet 

dimensionless streamline parameter 
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t 

u 

VXj Vy, Vz 


unit of tims^§ ^ , seconds 


velocity con 5 )onents of air, relative to free— stream 
velocity 


free— stream velocity of air, feet per second 

velocity components of a water drop, relative to free— 
stream velocity 


V 


w 


X, y, z 

y^iax 

a 

P 

A 

e 

0,r 


Xs 

d 

P 

Pi 


. Average velocity of a drop over an Interval of time, 
feet per second 

liquid water content of the air, grama per cubic meter 
«( 62 . 1 ^X 10 “®, lb /cu ft) 

rectangular coordinates, feet 

maximum ordinate of airfoil 

condition for which no drops will strike a body 

proportionality factor 

increment 

thickness parameter for transforming a circle into a 
Joukowski airfoil 

coordinates of .a circle of which the origin is at the 
center 

range of a drop in still air, feet 

absolute viscosity of air, pounds mass per second, foot 
density of a water drop, pounds per cubic foot 
density of air, pounds per cubic foot 
stream function 


Subscripts 

o depots s starting conditions where time is counted as zero 



MCA TN No. 1397 


5 


X 

denotes 

cr 

denotes 

n 

denotes 

n+i 

denotes 


one interval of tine 
a critical condition 
nth interval of time 
interval of time of (n+i) 


ANALYSIS 

To obtain area and concentration of water-drop impingement at 
the leading edge of ein airfoil, the path of the drops as affected 
by the air stream must be traced. In order to compute these paths, 
or trajectories, it is necessary to define the motion of a drop 
through the air. The equations which depict this motion, and which 
have been modified to include the effect of a drop depen*ting from 
Stokes* law of resisteince have been shown by Langmuir and Blodgett in 
reference 3 "to take the following form: 

(vi-u:c) (1) 


(Vy-Uy) (2) 


(^) = (3) 


dx 


2k 


Kv. 


dv- 


y = 


dy 


Ci^ 

2k 


dx 

= — Vy 

dt 


dy._ ■ 
d^ ■ ""y 


X T3 




The dimensionless quantity K, 
is evaluated by the equation 


used by Glauert (reference 2), 



w 



6 


NACA. TN No. 139? 


where in a system of consistent units^ p is density of the drop^ 
Pi is density of air, ft is radius of the drop, U is free— stream 
velocity, p is absolute viscosity of air, and c is the chord of 
the airfoil. 

It can he seen from equations (l) and (2) that 


dt 


2l^K 


(vx-ux) 


dt 


CpR 


(Vy-Uy) 


Since the average velocity of a drop over a s ma ll interval of time 
may he expressed hy 


~ _ Vn+idVn 
2 


and the average acceleration hy 


a = ^n+i'^n 
2 


the velocity and position of a drop at the end of any interval of 
time may he expressed hy the following equations; 


^^(n+i) 

(n+i) 




"'yn 


,24k ^ ^ j 

— /v^ — Uy N At 

24k yn ^n j 


(5) 


( 6 ) 
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*(n+i) = Vxn At - i 


~ (vxn-uxn) j + xn (?) 


y(n+i) = vyn At - I [ ^ (vyn-'^yn) 1 + yn 


(8) 


where the subscript n denotes a point at the beginning of some 
interval of time, and the subscript n+i denotes a point at the 
end of that interval of time. Ihe main assxunption made in these 
equations is that the acceleration of the drop is constant over the 
interval of time selected. 

Equations (3), (5)^ (6), (7), and (8) may be used to calculate 
the trajectory of a drop by a step-by-step process, provided that 
a sufficient range in values of air— velocity components ux, Uy, arxi 
the . drag coefficient function CI)R/2l^ are available. When Stokes* 
law of resistance holds, the value of Qi)E/ 24 is unity. The 
quantities Ux, uy, and C])R/24, are best used in large-scale plots 
to insure a fair degree of accuracy in the computations. In addition 
to the large-scale plots required, certain assumptions must be made 
in regard to the conditions of a drop at the start of its trajectory. 
This information is needed in order to continue with the step-by— step 
computations. It appears that the simplest manner in which to start 
the computations is by assuming the velocity of the drop to be the 
same as that of the air at a fairly large distance forward of the 
airfoil, and then estimating the velocity differentials Vx-ux and 
■'^y^y between the air and the drop at a new point. An estimation of 
the magnitude of these differences may be made by the following 
equations which are derived in appendix A: 


(uxq+®'^) i / {uxo+BAx) 2— 4Bux,Ax 

j ^ (9) 


^ -(BiV-Uyg) ± (Bily-UyQ)“-.4BUy^^ 
^yi " 


2 


( 10 ) 
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If 7o is the ordinate of a trajectory for a given drop size at 
infinity, the quantity of vater per unit time and area which progresses 

toward a hody is 

UwLyo 

I 

Since this aamfi amount of water distributes itself over the 
leading edge of the body, then. 


UwLyQ. = MLs 


and for a small increment 


M = Uw 


ds 


( 11 ) 


Since dyo/ds is proportional to the concentration of water impinge- 
ment at various points on a body in a water cloud, it can be appro- 
priately called a concentration factor. Thus, , , 

1 . . 

M = UwC (12) 


The total water collected by a body is the summation of all the 
degrees of concentration of water-drop Impingement over the leading 
edge of a body for a unit length of span such that 


Mp = f M|s = Uwy ^a-ic E (13) 

•Jo 

In equation (I3), E .is efficiency of water disposition over the 
leading edge, and is defined in terms of the maximum ordinate of the 
body, so that 


^^limit 


(14) 


where Vn-. ^ is the ordinate of a starting trajectory which just 
impinges upon the airfoil, and is the maximum ordinate of 

l^e body. 


When drops are quite small, they have insufficient momentum to 
carry them to the- body, hence, there is a critical size of drop which 
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possesses just s\afficient energy to Impinge upon the body. This 
critical condition is given hy 


a = 




dx 


1 


^( 15 ) 


To d.etermine the critical . drop size (i.e., the critical K value), 
it is only necessary to evalviate the quantity du^/dx at the 
stagnation point of the airfoil. Because the quantity du^^/dx is 
a property of the streamlines aroimd a body at the stagnation point, 
it has different values dependent upon the body chosen. 


CGmJTATION PKOCEDUEE 

The example given herein to explain the computation procedure 
is for a symmetrical Jouhowski airfoil, 12 percent thick, at 0° angle 
of attack. This airfoil was chosen because the streamline functions 
coiild be calculated with relative ease, and because its profile 
closely simulated that of an NAGA 0012 airfoil. Trajectory data 
for an NACA 0012 section were needed in connection with thermal ice- 
prevention heat— transfer computations. 

Streamlines were obtained about the Joukowski airfoil for a 
region estimated to cover the paths of drops of particular interest. 
The equations used to define the contour of the streamlines and the 
values of velocity components along streamlines were derived from 
Glauert (reference 4). These equations are given in appendix C. 

Streamlines calculated for the Joukowski airfoil are shown in 
figure 1. In addition, the contour of the NACA 0012 airfoil is 
presented to show its comparison with the Joukowski profile over the 
region where it was expected that trajectories would hit for drop 
sizes typically found in nature. 

Air— velocity-component data were used in the forms shown in 
figures 2(a), 2(b), 2(c), 2(d), and 3. Figure 2(a) gives the 
X— component of air velocity in terms of the x— coordinate in the 
vicinity of the airfoil, while figure 2(b) presents the same parameters 
several chord lengths forward of the airfoil. Figure 2(c) depicts the 
y— component of velocity in terms of the y— coordinate in a region near 
the airfoil. However, because of" their small magnitudes, y— components 

^Derived in appendix B. 
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of velocity well forweird of the airfoil had to "be treated somewhat 
differently than x— velocity components in order. to obtain values 
suffipiently accurate. It was fovind that the y— velocity components 
showed a good degree of change with respect to a change in ratios of 
the coordinates at any point. It was therefore foimd desirable to 
plot the y— velocity components in the meinner shown in figures 2(d) and 
3. Both X— and y— velocity components were plotted so that their valvies 
could be obtained to four significant figures. 

Because a water-drop trajectory can cross streamlines, and 
because air-velocity components are conveniently plotted in terms of 
streamlines, it is inqierative to be able to interpolate between the 
streamlines to ascertain on what streamline a drop is located at a 
given time. InteiTpolation between streamlines can be accomplished by 
plotting points on streamlines in terms of ratios of the point 
coordinates. Figure ^ shows the relationship obtained when this was 
done. The simplicity of this method of interpolation lay in the fact 
that the streamlines displayed a logarithmic separation, which permits 
the interpolation to be carried on by laying a logarithmic scale 
properly across the streamlines. It should be observed that, as a 
trajectory approaches the ordinate axis (taken in the computations 
to originate at the stagnation point), the ratio of the coordinates 
asymptotically approaches infinity. Therefore, in performing the 
step— by— step computations, it is advisable to abandon this interpola- 
tion procedure when values of the coor^dinate ratio are of the order 
of 1200 percent. For small values of x/c forward of the airfoil, the 
value of the stream function can be obtained for the drop with the aid 
of the information given in figure 1. 

Values of the drop drag-coefficient function Ci)R/2k and E 
were obtained' from reference 3> plotted to a suitable scale for 
use in the computations. These values are given in table I. As the 
step— by— step process progresses, new values of R and Ci^/2k are • 
selected from the plots. ■ 

From equations (3), (5), ( 6 ), (?), and (8) it is apparent that 
trajectories can be calculated for both ranges of K and Ry values, 
where Rxj is defined by 


Eu = 


2aUpi 


However, because of limited time available for performing the 
calculations, a set of trajectories (corresponding to several K 
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values) was computed for only, an Eu value of 95*65- The value of 
% was chosen because it corresponded to certain airspeed, altitude, 
and drop— size conditions of particular interest. 

Starting coordinates of trajectories were arbitrarily chosen for 
selected K values. As an example, by means of equations (9) and (lO) 
for a K value of O.O58I and starting coordinates Xo =210 percent 
chord, yo =0.750 percent chord, it x-^as determined that 




- Vj.^ = 0.9979 


Vy^ = 1.5X10- = 


CpE 

2k 


= 1.001 


Eu = 95.65 


A sample calculation for obtaining these values is shown in 
appendix D. With these values as starting values, a continuation 
of the computations is shown in table II for the trajectory. This 
table expresses equations (3), (5), (6), (7), and (8) in a convenient 
conq)utational form and is used in conjunction v/^ith figures 1, 2(a), 
2(b), 2(c), 2(d), 3, and 4 . 


EESULTS 

From calculations to determine the trajectories of water drops, 
the area and distribution of water impingement over the leading 
edge of an airfoil can be obtained. In addition, the total amount of 
water per imit length of span can be calculated. 
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If trajectories for one Ey value are plotted in the manner 
shovn in figure the distance along the surface to the stagnation 
point included by a trajectory may be measured. Then by taking 
included distances from a large-scale plot similar to figure 5, a 
figure like that shown in figure 6 is obtained. The limiting trajec- 
tories. beyond which no drops will impinge upon the airfoil surface 
may be obtained by connecting the end points for the several computed 
trajectories. This is shoxm by a -broken line in figure 6. It is 
interesting to note in figure 6, that for the case of K=«o, the 
problem reduces to one of plotting the airfoil ordinates against the 
included siurface distances. 

By obtaining the slopes of the cvirves in figure 6, according to 
eqxiatibn (ll), the concentration factor C may be obtained at 
different points over the leading edge of the airfoil. This operation 
was performed graphically, and the resdlts are shown in figure 7* 

The ratio of the limit starting ordinates of trajectories 
yQlimit to the tnariffluTn ordinate of the airfoil defines the efficiency 
of drop impingement (equation, (l4)). Therefore, when such ratios are 
obtained by means of figure 6, and plotted with respect to their 
corresponding K values, an indication of the efficiency of water 
deposition on an airfoil can be obtained. The value of K corre- 
sponding to zero impingement was estimated to be approximately 
0.0009 by means of equation (15)* The value of dux/dx needed 
to estimate this value of K was determined graphically by using a 
large plot of figure 2(a). A presentation of the efficiency of 
water-drop impingement on the JoukOwski airfoil, for a value of 
Ru= 95*65, is shown in figure 8. 

As a matter of interest, the paths of several drops of different 
inertia but the same starting points were obtained in relation to 
the velocity field about the airfoil. These relations are shown in 
figures 9(a) and 9(b), and were obtained by noting from the computa- 
tions (as performed by the procedure suggested in table II) the- 
corresponding air— velocity components at different times during a 
drop trajectory . 


DISCUSSION 

The procedure demonstrated herein for calculating water-drop 
trajectories about a body in two-dimensional flow is convenient to 
use when aid from a differential anailyzer is not available. In the 
use of the method, however, certain errors are involved. One of the 
errors is involved in the ability to predict, with reasonable accuracy. 
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the starting values of the velocity differentials Vxj^— and 
Vyj^-Uyj^ for small drop sizes (K values less than approximately 
0.02). It appears that, in order to calculate trajectories for these 
small drops, either the air— velocity components must be obtained to 
great accuracy or else a better method of approximating the starting 
values of the differentials must be derived. Another error lies in 
the accumiilation, during successive integration steps, of errors 
Involved by the assumption of an avera^ constant acceleration over 
the time interval selected. Near the airfoil, errors accumulate more 
rapidly than they do well forward of the airfoil because of the 
increased rate at which values of the air— velocity components change. 
(See figs 9(a) and 9 (b).) Large changes in air— velocity components 
mean large changes in selected values of CpR/24 which multiply the 
velocity differentials v^-Ux and vy— Uy. Therefore, it is desirable 
to reduce the magnitude . of the arbitrary time interval when a 
trajectory approaches the leading edge. To determine a suitable time 
interval to apply in the computation of a trajectory as it approaches 
a body, a criterion is to permit the drop to travel, in the time 
interval selected, a distance in which the air— velocity x-con5)onent 
values do not change by more than approximately 10 percent. 

To greatly reduce the amount of work involved in con5>uting a 
set of trajectories, starting ordinates of trajectories can be taken 
very close to the body, eind time intervals can be chosen as large 
as is commensurate with the criterion for their selection. When this 
was done for the computation given in table II (the starting ordinate 
being taiken at half a chord length forward of the airfoil), area of 
water-drop impingement was reduced 15 percent. To start the calcula- 
tions so near the airfoil, it was assumed the drop followed the 
streamline up to the starting point. Because of this assumption the 
starting ordinate was large by approximately 10 percent, compared to 
2 percent for the computations given in table II. 

With the set of trajectories calculated by the step— by— step 
process and presented in figiare 6, a comparison can be made between 
the impingement of water drops on the Joukqwski airfoil and its 
ecLuivalent cylinder, for a small range in drop sizes. Such a comparison 
is made in figure 10, where the airfoil was used as the basis for 
con5)arison. For area of impingement, the equivalent cylinder can be 
used without much error to a drop— size diameter of approximately 20 
microns. However, for rate of impingement, appreciable difference 
is not encountered until drop diameters are above 3.^microns. It is 
of interest to note that meteorological data obtained in flight 
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(reference 5) indicate a mean-effective drop-size diameter^ for 
cumulus clouds of I5 to 25 microns. Furthermore, mean-effective drop 
diameters as high as 50 microns have been measured in clouds of the 
altostratus category. Since, in accordance with the definition of 
mean-effective drop diameter, half the water contained in a cloud 
sample is composed of drops larger than the mean-effective diameter, 
it is apparent that the use of an equivalent cylinder may he subject 
to considerable error when applied to typical icing conditions in 
which large drop— size diameters have been observed. 

The method for computing trajectories of water drops in two- - 
dimensional flow can be extended to the case of three-dimensional 
flow about a body of any shape. To extend the method, the trajectory 
of a drop must be expressed for the third plane. The additional 
trajectory equations needed are: 




dz 

— = 
dt 

In addition, equation (3) which connects the variable drag coefficient 
Cd with the local Beynolds number of the drop at a given time becomes: 


+ (vj.-Uj)^ (17) 


Thus when equation (I6) is placed in forms similar to equations ( 5 ), 
( 6 ), (T), and ( 8 ) in the analysis, and equation (I7) is solved in a 

^ Jifean-effective drop diameter as used in the present report is the 
size of drop in a cloud sample for which the amount of liquid 'water 
existing in water drops larger than that drop is equal, to the amount 
of liquid water existing in drops smaller than the drop. 
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manner eimilar to that indicated in table II, a three-dimensional case 
may be handled by step-by— step integration when air— velocity cpmponents 
and starting conditions are known. Application of trajectory computa- 
tions for three-dimensional flow can be made to the case of windshield 
ice prevention, where it is desirable to know the rate at which water 
impinges upon a windshield in determining heat requirements to keep 
a windshield clear of ice. Equations (I 6 ) and (17) are basic ones 
which can be adapted also for use with a differential analyzer. 

The, concentration factor C in equation (12) for the three- 
dimensional case is now evaluated by . . 


dAp 

dAs 


( 18 ) 


In equation (I 8 ), Aq is an area which is outlined by several 
starting trajectories; and Ag is the area, on the surface of the 
body, which is outlined when these trajectories intercept the body. 

The rate of water impingement My in equation (12) is evaluated 
for the three-dimensional case by 


Mr = UiiffiAfrontal ( 19 ) 


In equation (19) > E, the efficiency of water-drop deposition on the 
body, is now expressed by • ' 


,jj _ ■^oiimit 

Afrontal . 


where is the area included by starting trajectories which 

just impinge upon the body, and Afrontal is the frontal area of the 
body. 
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As has already been mentioned, the limit distribution of vater— 
drop impingement on the leading edge of a body can be obtained from 
the body ordinates. It is interesting to observe that, for yery high 
speeds, the limit condition K=oo is approached. Hence, it is 
believed that a comparison of the water-drop impingement over several 
different airfoils for the limit condition is of interest. Such a 
con5>arison is made in figure 11. For these limit conditions, any body 
having a leading— edge radius will have a value of the concentration 
factor at the leading edge of unity; whereas, any body without a 
leading-edge radius will have a much lower value of the concentration 
factor depending upon the slope of the contour near the stagnation 
point. According to equation (I3), the area under the curves in 
figure 11 is proportional to rate of icing per unit span, and in 
accord with theory, for airfoils with the seimc 'value of maximum 
thickness^ regardless of location, the areas under the curves in 
figure 11 -are equal. The difference between each of the three 
airfoils lies in the manner in which water drops distribute themselves 
over the airfoil. The reason for this difference is because a limit- 
ing water— impingement distribution corresponds to the case where drops 
are either so large or,. are traveling so fast that they are undeflected 
by forces in the air stream. Therefore, the effect of shifting the 
location of maximum thickness is to change the degree and position of 
water-drop impingement. Thus the limit water— drop— impingement distri- 
bution for a given body can be made to have any desired profile by the 
proper combination of contour and location of maximimi thickness. 

CONCLUSIONS 

From the method outlined in this report to obtain trajectories 
of water drops forward of a body, the following conclusions can be 
drawn; 


1. By a step-by-step process which considers deviations from 
Stokes* law of re si stance, the area and distribution of water-drop 
impingement on a body of any shape and angle of attack can be calculated 
when the air— velocity components are known. 

2. The rate and area of water-drop impingement on an NACA 0012 
airfoil having an 8-foot chord, is approximated with reasonable 
accviracy by the impingement on a cylinder of radius equal to the 
leading-edge radius of the airfoil for only drop diameters below 20 
microns at speeds of about 200 miles per hour and at an altitude of 
7000 feet. 
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3. The critical drop— size diameter corresponding to zero impinge 
ment can be determined for a two-dimensional body if the velocity 
components are kiiown at the stagnation point. 

4. The area and distribution of water-drop impingement on 'an 
airfoil in motion through a cloud of known drop— size distribution can 
be obtained by calculating the area and distribution of impingement 
for each drop— size range present. 

. 5* Thin airfoils of different sections (such as wedge— shape^ 

NACA 0006, and circular— arc airfoils), having the same exposed frontal 
areas, will have approximately the sane rates of water-drop impingement 
at hi^ speeds, but the distribution of water-drop impingement will 
be different for each section, and will- depend upon the rate of change 
of the slope of the airfoil surface and the location of wayi mum 
thickness. • 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. ^ July 194 7 
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APPENDIX A 

Derivation of Starting-Point Equations 

. At a position Xot yo well forward of the airfoil, it is 
assumed that a drop is moving with a free— stream velocity u^. Then 


UXo = ^Xo (Al) 

^yo = "^yo 


The position of the drop after an interval of time At is 


Xx = Xo + Ax 


yi = yo + ^ 


where Ay may he expressed in terms of Ax hy 


Ay = 


^yp 

UXo 


(Ax) 


(A2) 


The velocity of a drop, starting initially from a point Xq, y^, 
after an interval of time At is 


vxi = VXq + Avx 


(A3) 


''yi - ’'yo * '^’'y 
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From tasic trajectory equations 


Vx <ivx = B (vx-ux) <ix 


(Ah) 


where 


B = 


QpR 

2 UK 


By substituting equations (Al) and (A3) into equation (a 4 ) 


Vxi (vxi-%o^ = ® (vxi-«xi) ^ ■ 


and by solving for vx^. 


vxx = 


(u^^+BAx) ± \/ (u3j^^+BAx)^-4Bu3j.^ Ax 

I, I ,1 

2 


(A 5 ) 


In a manner similar to that by which equation (A 5 ) was derived, 
the initial value of the y-component of the drop velocity v„ may be 
obtained. The value of Vyj^ becomes 


vyi =' 


- (BAy-Uy^j) (BAy-Uy^ )2n4Buy^Ay 


(A 6 ) 


2 
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APEENDIX B 

Derivation of Equations for the Condition Where 
No Drops Inq)inge Upon a Body 

In deriving the relationships which describe the limiting value 
of E below which no drops will impinge upon a body, it is necessary 
to consider the basic trajectory equations and the flow field at the 
stagnation point of the body. 

Near the stagnation point of a body, the velocity of both air and 
of drops which nearly follow the :streamlines is so small that CpR/24 
may be taken as unity. In addition, the air— velocity components are 
nearly linear in relation to their distance from the stagnation point. 
Glauert shows in reference 2 that K is proportional to the range of 
a drop in still air; and in reference 3 } Langmuir denotes this range 
by Xs, where 

K = ^ 
c 


Since K is a function of the air— velocity components and 
because^ near the stagnation point of the body x— components of air 
velocity according to reference 3» simplified to 


ux » ap 


(Bl) 


the range of a drop may be written in terms of a distance 
p forward of the stagnation point of a bodj’^. This distance is 
defined by 


P = 


cx 

>^s 


Equation (l) may be written 


dvx 

Kvx = Vxr<iP 


(B2) 


(B3) 



MCA TN No. 1397 


21 


To solve equation (B 3 ), v^ may "be replaced in manner similar 
to ux "by 

Vx = 3 p (B4) 


where by differentiating 


dvx = Pdp + pd3 


Then 


= 3 +P-2. 

dp dp 


(B5) 


By differentiating equation (B2) with respect to x. 


Kdp = dx 


(B 6 ) 


by substituting equations (BU) and (B6) into equation (B 3 ), 



(B7) 


By setting equations (B5) and (B 7 ) equal to each other and rearranging 
terms, 

^ = _ (b 8 ) 

p Or-3+3^ ' 


For the integration of equation (b 8 ), Pierce’s tables (reference 6 ) 

gives 


p p(ip ^ 

.] ct-3+3 2h 2h J X 


(b9) 
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where 


X - a+g3+h3^ 


and 



tan”^ 


f agtgA 

vyr J 


In eq^uation (B9) 


s Uahr-g2 


Since g=-l and h=l, equation (B8) Becomes 


Inp - -ln(or-3+P^) r- .. tan i ^ 

^ -v/Uor-l VVUc^l/ 


(BIO) 


In equation. (BIO), it is apparent that a= — is a critical value 
in evaluating the function of p. 

A substitution of equation (Bl) into equation (B2) gives 


Ku^ = cue 


a nd by differentiating with respect to x,. 


Kcr 



a 


Since equation (BIO) is dependent 'solely upon the trajectory 
equations, the value of a may always be taken as l/U in evaluating 
the critical K value below which no drops will inq>inge on a body. 
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Eo^jovoi', the quantitj dux/dx is dependent upon the shape of a body, 
the attitude of the body in the air stream, euid the type of flow around 
the body. It, therefore, must be evaluated at the stagnation point for 
the flow around the body under consideration. 


APPENDIX C 

Equations for Obtaining Streamlines and Velocity 
Components about a Symmetrical Joukowski 
Airfoil, 12 Percent Thick 

The contour of a streamline is given by 


X = (r cos d + c) ( 1 + ^ ) 

r^+€?+2re cos 0 / 


. y = (r sin 0)( 1 : r — i ) 

V r^4e^+2re cos 0 / 


in which 


r 


s 


2 sin 0 


V \2 sin 0/ 


(Cl) 


In equation (Cil) 



where \(f , a stream function, can have any desired value. In addition, 
for a symmetrical Joukowski airfoil, 12 percent thick, the value of e 
is 0.1021. 
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The components of air velocity in terms of free— stream velocity 
are given hy 


= U1U2 — V1V2 


where 


Uy = U1V2 + usvx . 


. I 

^ cos 29 

-r^ ■ 

(1+e)^ 

= sin 2a 


ua . r 1 

L (D-l)^+F^ J 


(C2) 


v^ 


F 

(D-1)2+F2 


(C3) 


Inequations (C2) arid (C3), 


D = cos 20 + 2 r€ cos 0 + e® 


F = 2r sin 0 (r cos 0 +'€ ) 
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APPENDIX D 

Computation of Starting Values for the 
Step-by-Step Computations. 

According to equation (A5), the x— velocity component of a drop is 


vxi 


_( — 4B Ux 2 .^^ 


In order to obtain starting values with this equation, this equation 
must be solved by a trial— and— error method. A procedure for a trial— 
and-error solution is demonstrated for the starting values used in 
table II . 

For the conditions of table ll. 


Xo = 2.1; yo = 0.0075 


K = O.O58I; Eu' = 95-65 


Several chord lengths forward of the airfoil, as a first approx- 
imation, the quantity vy^^— uyj^ may be considered zero. Thus 


R«±Eu (vxi-uxi) 

Assuming a value of— (vx^— Ux^) = 0.0001, the value of R is 

R - 0.00956 
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for which the value of CpR/2h is 


CdR 

« 1.001 

oil 


With Xo as 2.1 chord lengths forward of the airfoil, and Xi as 
2.0, the corresponding magnitudes of air— velocity components are 
ux = 0.9979 and uxi=0.9978. The value of Ax is then 0.1,‘and the 
va£ue of B is 17 . 228 , from which the value of vx^ from equation 
(A5) is 0 . 9979 . Therefore, the initial assuuQition of 
-(vxi“iixi) =0*0001 was very nearly correct. 

Since the quantity (vyj— Uy^^) contributes very little to the 

initial value of CjjR/ 24, once the approximate value of Vx^-Ux^ 
has been determined, an approximation of Vy^-Uyj^ may be made by 
equation (A6) without changing the initial value chosen for CdR/ 24. 
For the trajectory calculated in table II, yyi“Uyi was assumed to be 
0.1xl0~®, and it was also found by equation (A2) and (A6) to be of 
this magnitude. 

After the differentials ^yi“^i evaluated, 

they are entered in the appropriate columns in table II, and the 
process continued by the steps indicated in the table. 
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lABLB II— COMPUTATION OP A W4TEE-DR0P TRAJECTORY 
[RU® = 9,15 X 10®| K = 0.0681] 
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TAa.B II.- COBTnroKD. COamATIOT O? A miSR-IfflOP THAJKTOBT 
- 9.16 X 10 »t K . 0.06811 
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T*Ht.R II,- CONTIITOED, COMPOTATION OP A WATER-DROP TRAJBCTQRT 
I Hd* » 9.16 X 10*1 K ■ 0.0681 ] 
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PERCENT THICK SHOWING A COMPARISON BETWEEN THE SURFACE 
5TREAMUNE5 FOR AN MCA 0012 AND THE JOUKOW5KI AIRFOIL 
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Fig. 2c 
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NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 

FIGURE 2C.- AIR-VEIOCITY-COMPONENT PROFILES FOR VARIOUS 
5TREAHL1NES ABOUT A SYMMETRICAL JOUKOW5KI 
AIRFOIL (2 PERCENT THICK. Y' VELOCITY COMPONENTS 
IN PROXIMITY TO THE AIRFOIL 
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Fig. 2d 
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1.2 M 6 a >0 20 MO 60 80 lOO 200 

RATIO OF Y-COOROINATE TO X- COORDINATE, IN PERCENT 

FIGURE 2D.- AIR-VELOCITY COMPONENT PROFILES FOR V/«10U5 STREAMUNE5 ABOUT A SYMMETRICAL 
JOUKOWSKl AIRFOIL 12 PERCENT THICK. Y- VELOCITY COMPONENTS IN TERMS OF 
COORDINATE RATIOS 
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Fig. 3 



FIGUHC 3.- COORDINATE RATIOS OF POINTS ON STREAMLINES AS A FUNaiON OF Y- COMPONENTS 
OF AIR VELOCITY FOR. A. SYMMETRICAL JOUKOWSKI AIRFOIL 12 PERCENT THICK. 
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distance along chord datum L\NL \N PERCENT 

FIGURE H.- STREAMUNtS ABOUT A SYMMETR\CAL J0UK0W5K1 AIRFOIL 12 PERCENT THICK 
AS A EUNTCTION OE THEIR COORDINATES. 
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Fig. 4 cone. 
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FIGURE 5.- RELATION OF WATER-DROP TRAJECTORIES TO STREAMLINE 
FIELD ABOUT A SYMMETRICAL J0UK0W5K1 AIRFOIL 12 
PERCENT THICK. 
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Fig. 6 
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DISTANCE ALONG SURFACE OF AIEPDIL FROM LEADING EDGE, ^ IN PERCENT 
FIGURE 6.- AREA OF DROP INTERCEPTION CORRESPONDING TO DROPS OF DIFFERENT INERTIA 
FDR A SYMMETRICAL JOUKOWSK.I AIRFOIL IZ PERCENT THICK.. 
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Fig. 8 
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EFFIQENCV OF DROP IMPINGEMENT, E, IN PERCENT 
FIGURE 8 .- WATER -DROP -COLLECTION EFFICIENCV FOR A SYMMETRICAL JOUlCOU/SKI AIRFOIL 
(2 PERCENT THICVC. 
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Fig. 9b 



DISTANCE FROM CHORD DATUM L1NE,Y^^ PERCENT 

NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 

FIGURE C|B.- influence T)F VELOCITY FIELD ON THE PATH OF WATER 
DROPS WMm DIFFERENT INERTIA FOR A SYMMETRICAL 
JOU<OW5Kl AIRFOIL 12 PERCENT THICK. Y- VELOCITY 
COMPONENTS 
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